arXiv: 1501.07023v2 [hep-ph] 4 Mar 2015 


On the nature of ai(1420) 

M. Mikhasenko, 1 B. Ketzer, 1 and A. Sarantsev 1,2 

1 Universitat Bonn, Helmholtz-Institut fur Strahlen- und Kernphysik, 53115 Bonn, Germany 
2 Petersburg Nuclear Physics Institute, Gatchina, Russia 

The resonance-like signal with axial-vector quantum numbers J pc = 1 ++ at a 
mass of 1420 MeV and a width of 140 MeV, recently observed by the COMPASS 
and VES experiments in the /o(980)7t final state and tentatively called ai(1420), is 
discussed. Instead of a genuine new meson, we interpret this signal as a dynamical 
effect due to a singularity (branching point) in the triangle diagram formed by the 
processes ai(1260) -* I\*K, K* —> Kn. and KK /o(980) (+ c.c). The amplitude 
for this diagram is calculated. The result exhibits a peak in the intensity with a 
sharp phase motion with respect to the dominant ai(1260) —> pir S- wave decay, in 
good agreement with the data. The branching ratio of ai(1260) —»• /o(980)7r via the 
triangle diagram is estimated and compared to the dominant decay ai(1260) — > pir. 
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I. INTRODUCTION 

Understanding the quark interaction at low and intermediate energies is one of the most 
challenging tasks of the theory of the strong interactions. To create a theoretical non- 
perturbative approach or, at least, to build a reliable model one should understand the 
nature of strongly interacting particles and their excitation spectrum. The classical quark 
model assumes that mesons are bound states of quarks and antiquarks, and groups the 
low-mass states into nonets with the same spin J, parity P and charge-conjugation parity 
C, with fixed mass differences between nonet members. The second assumption is that 
the quark-antiquark interaction at large distances is governed by a linearly rising potential 
which explains the phenomenon of quark confinement and predicts the full spectrum of 
quark-antiquark excited states. The success of the quark model is indisputable: most of the 
known mesons correspond very well to the predicted scheme (TJ. 

However, it seems that the meson spectrum is notably richer than that predicted by the 
quark model. There is a growing set of experimental observations of resonance-like structures 
in partial waves with quantum numbers which are forbidden for the quark-antiquark system 
or situated at masses which can not be explained by the quark-antiquark model, see e.g. 
[21 El and references therein. 

Recently the COMPASS [1H6] experiment reported the observation of a small resonance¬ 
like signal with axial-vector quantum numbers I G (J PC ) = 1“(1 ++ ) in the /o( 980 ) 7 r P -wave 
of the 7r _ 7r _ 7r + final state, produced by diffractive scattering of a 190 GeV n~ beam on a 
proton target. The signal was also confirmed by the VES experiment [7] in the 7r _ 7r°7r° final 
state. In both experiments, the three-pion final states were analyzed using a two-step partial- 
wave analysis (PWA) technique. In the first step the data were grouped in small bins of 
37r-invariant mass and momentum transfer. The isobar model was employed to parametrize 
possible decays to three final pions. An isolated, relatively narrow peak was found in the 
intensity of the 1 _ 1 ++ / 0 7rP-wave at a mass around 1.4 GeV, accompanied by a sharp phase 
motion of this wave relative to other known resonances, with a phase variation exceeding 
180°. In the second step of the COMPASS and VES analyses the spin-density matrix 
resulting from the first step was fitted with a model including Breit-Wigner resonances 
and background contributions. The new signal was described rather well with a hitherto 
unknown resonance, which was tentatively called ai(1420) with a mass M ai ~ 1.42 GeV and 
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width r oi ~ 0.14GeV. The interpretation of this signal as a new state in the framework of 
the quark model is difficult, ft cannot be a radial excitation of ai(1260) which is expected 
to have a mass above 1650 MeV. It is also not expected that the radial excitation has a 
width which is much smaller than the one of the ground state. Therefore, this signal is 
to be considered either as a strong candidate for a four-quark bound state or a meson- 
meson molecular bound state or to be explained as a some dynamical effect resulting from 
multi-particle interaction. 

In the present paper we show that a signal of comparable strength, including the rapid 
phase motion, can be expected by the opening of the K*K decay channel for the isovec¬ 
tor ai(1260), and the re-scattering of the kaons. There are two isospin combinations of 
intermediate particles: 

(i) ah (1260) -)■ K*°K~ -> vr ~K+K~ n~f 0 , 

(ii) ah(1260) -)■ K*~K° -> n~K 0 K° -> vr ~f 0 . 


The corresponding triangle diagram has a logarithmic singularity on the tail of the wide 
ai(1260) which is due to a very peculiar kinematic situation, in which all intermediate 
particles are almost on their mass shell, causing a resonance-like effect. 

Recently, Wu et al. showed that the same triangle singularity in the isospin-0 chan¬ 
nel can account for the anomalously large isospin violation effects observed by BESIII for 
?7(1407/1475) and /i(1420) —» /o(980)7t — y 3n [10]. A similar dynamic mechanism was sug¬ 
gested earlier by Achasow and Kozhevnikov m to explain the resonance-like signal observed 
in the cftn 0 mass spectrum of the reaction iT~p — y (ftir°n [T2] by the decay of p(1700) — y K*K 
and rescattering of KK — > (ft. Triangle singularities are currently also being discussed in the 
context of the newly discovered XYZ quarkonium peaks m- 

Our paper is organized as follows. In sections |TT| to m we only discuss the triangle dia¬ 
gram for process (i) with intermediate particles (K*°,K + ,K~), the calculation for process 
(ii) proceeds analogously. In section [H] the kinematic conditions for the appearance of the 
triangle singularity are analyzed. The amplitude for the triangle process is calculated in the 
following two sections. In section [HI] we first present an approach to calculate the imaginary 
part of the amplitude making use of Cutkosky cutting rules and the calculation of discon¬ 
tinuities. This method helps to understand the structure of the amplitude singularities. In 
section IV we then use an effective Lagrangian approach |1][J to calculate the full amplitude, 
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i.e. the real and imaginary parts, needed to predict the phase motion. In both sections, we 
start with the case of scalar particles to illustrate the underlying physics. In scalar theory 
the behavior of the amplitude is oc log(s — E \) near the singularity. Since the amplitude 
behavior could be different in scalar theory and in interactions of particles with spin mi, 
then the realistic situation for particles with spin and finite width is calculated. We show 
that the singularity is removed only by including the finite width of unstable particles and 
gives a contribution r\_/ log r k * ■ 

In section |v| we estimate the branching ratio of ai(1260) —* / 0 ( 980 ) 7 t _ , now including 
both isospin combinations, relative to its dominant decay channel pn~, and compare the 
signal we expect for the triangle diagram to experimental values as reported in [6] . 
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II. KINEMATIC CONDITIONS FOR TRIANGLE SINGULARITY 



Figure 1. The process 0 —> 1 + 2 for particles with 4-momenta po, Pi, P 2 , proceeding via a triangle 
diagram with intermediate particle momenta k\, k 2 , k^. 

It is well known that logarithmic singularities arise in processes which proceed via the 
triangle diagram shown in Fig. [l] As it was shown by a general analysis of singularities in 
scalar theory [TS], the amplitude behavior near the branching point of a cut is oc log(s —So), 
where s is an external invariant. The position of the singularity s 0 can be obtained from 
the simple condition that all intermediate particles are on mass shell and collincar to each 
other. It is given by the system of Lorentz-invariant Landau equations: 

' 

kf = m 2 , z = 1... 3, 

< x k Xfl -yk 2fl +zk 3fl = 0, x, y, z G [0,1], (1) 

x + y + z — 1, 

with ki and rrii the 4-momenta and masses of intermediate particles, respectively, and x, y, z 
the so-called Feynman parameters. The system of equations ([Tj) for x, y, z is overdetermined, 
so it is solvable only in exceptional cases. For the special case of the decay of a] - (1260) to 
/o(980)7t~ through intermediate particles (K*~,K°,K°) or (K*°,K + ,K~), and neglecting 
the finite width of the / 0 (980), the external momenta pi and p 2 depend only on s = p% (see 
Fig. 0 for the definition of pi). Using kinematical relations between internal and external 
momenta it can be shown that the system ([l]) has solutions only if yfs = Fd, 2 , where 
Ei = 1.42 GeV, E 2 = 1.46 GeV. These pinch singularities are shown as dots in Fig. [2} As 
can be seen, the conditions x,y,z G [0,1] are satisfied only for the first solution. 

Here we give a simple kinematic explanation for the appearance of the singularity. The 
initial state ai(1260) with J PC = 1 ++ can decay to real K*K + c.c. starting from the 
threshold y/s = 1.39 GeV. Then the K* decays to real K and tt. Note that the K from K* 
decay can go to the same direction as the K, the ratio of velocities of K and K is a function 
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Figure 2. Diagram illustrating the positions of singularities of the triangle diagram shown in Fig. [l] 
for the Feynman variables y, z. The light-gray triangle is the kinematically allowed region. The 
black dots are pinch singularities, corresponding to y/ s = E\ (lower dot) and y/s = E 2 (upper dot), 
respectively. The curves are solutions of A yz + rn\ (1 — y — z) = 0 (see equations 
y/s = E\ (solid line) and y/s = E 2 (dashed line). 


16 


and 


17) for 


of y/s as well as of the invariant mass of K and K. The invariant mass of K K going in the 
same direction is equal to the mass of / 0 only if y/s = 2 , but for E 2 the K is faster than 

K and thus the K cannot catch up the K to form / 0 . Only for the solution E] do the K 
and K proceed in the same direction with the same velocity with an invariant mass equal 
to that of the / 0 . 

The kinematics discussed here demonstrates a very peculiar situation in the decay of the 
ai(1260) to K*K + c.c: just above the two-body threshold, the re-scattering in the triangle 
can happen with particles almost on mass shell. 
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III. IMAGINARY PART OF THE AMPLITUDE 



vK+(k 3 ) 

/o(980) 


Pi 


P2 



Figure 3. Two possible cuts which contribute to the imaginary part of the matrix element of the 
process oi(1260) —> /o(980)7r“. 


In order to understand the structure of the amplitude, we first consider the imaginary part 
only, based on discontinuities. The technique was developed by Cutkosky [16j, is described 
e.g. in the Gribov lectures ra. and was successfully applied by Achasov and Kozhevnikov 
for similar process, p' —> (pn HU- The imaginary part of the amplitude M of the diagram in 
Fig.0 is related to the discontinuity across the cuts shown in Fig. [3] by 


ImM B1 ^ /of = - (Disc K * A - + Disc kR ) 

To calculate the discontinuities, we use the following expression 

d 3 h 


( 2 ) 


Disc = 


n 


(2ir) 3 2 E* 


X 


l J2 -Ms) 

\polarization / 


x (27r) 4 m(mom. cons.) , 


(3) 


cut v ' % \polarization 

where are matrix elements for processes on the left and right hand side of the cutting 
line, respectively (see Fig. [3]). We are calling particles which are crossed by cut line as cut 
particles. The integration is over all momentum space for cut particles, i.e. k j are momenta 
of cut particles, Ef are the corresponding energies. If a cut particle has spin we sum over 
all possible polarizations. 


A. Simple model with scalar intermediate particles 


For the case of scalar intermediate particles, the expressions for the discontinuities are: 

d 3 /ci d 3 /c 2 1 


Di sc ( k*k = 9 3 


(2vr) 3 2Ef (2vr) 3 2 E$ m 3 — + ie 


- x (27r) 4 (f 4 (po - h - k 2 ) , (4) 


x 















T^v * (sc) 3 

Disc// = 9 


d?k 2 


d :i b. 


x 


x (2 tt) 4 (T(fc 3 + k 2 -p 2 ) • (5) 


{2tt) 3 2E^ (2tt) 3 2£| m\-k\- ie 

Here, the products of matrix elements Mi • are given by the coupling constants at the 

three vertices, which are set to g , and the propagator, which is a function of the angle between 
k] and p\ in equation (Jd| and a function of the angle between k 2 and p\ in equation ([ 5 ]). 
For both discontinuities, the cut particles ( K *, K~ and K + , K~, respectively) are set on 
their mass shells. Integration with delta function in equation ([5]) is performed in the fo rest 
frame. After carrying out the integration we arrive at the following expression: 


Iml 


ir( sc ) 

“■<21—>/07T 


167T 


, A + l + ie 
log -- : -- + 


2 |p | \/s A-l + ie 2\p'\M-[ 


log 


C + 1 — ie 
C - 1 -ie 


where the coefficients A, C originate from the propagators, 


A = (ml 


ml-Ml + 2E\E 1 l)/(2\k\\p\) 


C = (ml - s - ml + 2E' 0 E?)/(2\k'\\p'\) 


( 6 ) 

(7) 

( 8 ) 


Here, Mf = p 2 , \k\ = / (2yfs), \p\ = A 1 / 2 (s, Mf, M|)/(2 v /s) are the momenta 

of the corresponding particles in Gq rest frame, with the Kallen function 


X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + yz + zx) 


(9) 


E\ = (s+ml — ml)/(2 t/s), E\ = (s+Mf— M|)/(2y / s). The values with prime are calculated 
at the fo rest frame: \k’\ = A 1/,2 (M|, m|, mf)/(2M 2 ), \p'\ = A 1 / 2 (s, Mf, M|)/(2M 2 ) are 
the momenta of K + (K~) and 7r _ (ab)- The corresponding expressions for energy are 
E% = (Mf + ml — m|)/(2M 2 ), E' 0 = (s + Mf — Mf)/(2M 2 ). The imaginary parts of the 
expressions Q. ([5]) compensate each other and the equation ([6]) is real. 

The imaginary part of the amplitude ImM^ i /j.j 0 ^('S) as well as the contributions from 
the individual discontinuities are shown in Fig. [4] by dashed lines. One can clearly see 
singularity a/s = E\ and E 2 values. The singularity in E 2 is out of kinematically allowed 
region of reaction, so the sum of the two discontinuities is smooth at E 2 . One can also 
notice that the imaginary part is not zero below the threshold of K*K threshold. Here, 
the contribution comes from Disc A -^, because the mass rrif 0 is above 2 m K ± threshold. Of 
course, taking into account the real shape of K* and fo will make the amplitude smoother, 


as shown in section IV C but the effect of the singularity at a/s = E\ will remain. This 


conclusion will also not change when the spin of the particles is taken into account, as will 
be shown in the next section. 
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B. Realistic case: VPP intermediate particles 


In reality, the particles involved in the process carry quantum numbers different from the 
scalar particles used in the previous section. The ai(1260) with axial-vector quantum num¬ 
bers J p = 1 + decays to K*K with vector and pseudoscalar quantum numbers, respectively. 
The K* decays to two pseudoscalars, Kn. The Feynman rules for the hadronic vertices 
which we use are given in appendix [A} 

The expressions for the discontinuities are: 


Disc 


(vpp) 

k*k y 


3 


Disc 


(vp_P) 

KK 



d 3 Aq d 3 k 2 

(2'Kf2E\ (2vr) 3 2 E\ 

d 3 k 2 d 3 k 3 

(2vr) 3 2 E\ (2vr) 3 2£ 3 fc 


( g pv - (pi - h) 

m 3 — k'l + ie 

£ov (.9^ ~ (Pi - h) 
ml — k\ — ie 


x (27r) 4 5 4 (po - ki - k 2 ), 

( 10 ) 

x (27r) 4 <5 4 (fc 3 + k 2 -p 2 )- 

( 11 ) 


Here Eq is a polarization vector of the cq state, we use notation shown in Fig. [3j 


Im, Disc^ K l 2, Disc*- K l2 



Figure 4. Energy dependence of Im M^ c 2^ o7r (s) and ImA^j oJr (s) (dashed and full lines, respec¬ 
tively). The contributions of the discontinuities K*K and KK are shown by green and blue lines, 
respectively. 


After integration, we have: 


Disc 


(vpp) 

k*k y 


3 


1 2\k\ 

8tt t/s 


(eo(pi ~p 2 ))x 


Mf + m\ — ml 
Arrifpl 


1 + 


( M m l , A-l + ie 

\ \k\Ml + m\-ml 2 J g A + l + ie 


x 


( 12 ) 
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Disc 


(vp_p) 

KK 


= ~9 


1 2\k'\ 


x 


87T M2 

M{ + ml — ml E' 0 


Mpi ~P2))x 


4mfp 


12 


Mo 


1 + 


\p'\ 2m\M 2 — (Mf + rn\ — m^E^' C 


\k'\ 


(.Mf + ml — ml)E' 0 


+ 7T n°g 


C-l-ie 


C+l 


re 


(13) 


The notations could be found in the previous section. The P -wave from K* decay is prop¬ 
agated to the /o7T P- wave. So we have a factor (e 0 (pi ~ P 2 )) in the final expression for the 
imaginary part of the matrix element. We separate it to compare the result with the scalar 
case: 

= 9 3 ^^(eoiPi ~ P 2 )) ■ (14) 

A 0)^/0 ?r i s plotted in Fig. [d] together with the result from the scalar theory. The two results 
are very similar. 
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IV. FULL AMPLITUDE FOR a* (1260) -> / 0 (980) ir~ VIA K*°K + K~ TRIANGLE 


After the calculation of the imaginary part of the amplitude based on discontinuities we 
proceed now to the calculation of the full amplitude for the triangle diagram shown in Fig. [5] 
using Feynman rules for hadronic processes in an effective Lagrangian approach ia (see 
appendix [A] for the parameterization of vertices). As in the previous section, we start from 


the simple case of scalar particles, and generalize to particles with spin in section IV B 


Po 


K*°(ki) 
cq(1260) 


K~(k 2 ) 



/o(980) 


Figure 5. ai(1320) —> fo(980)ir triangle diagram. 


A. Simple model with scalar intermediate particles 


In case of vertices involving scalar particles only the matrix element for the triangle 
diagram in Fig. [5] is 

d 4 £q 1 


I (sc) - a 3 
V-UoTT - y 


(2tt)H (mf — kf — ie)(m| — (p 0 — ki) 2 — ie)(m| — (hi — p ±) 2 — ie) ’ 

We calculate the integral using the standard technique of Feynman parameters and Wick 
rotation, which is described in more detail in appendix |B| 

1 


ir(sc) _ 9 


fl-y 


167T 2 


d y 


d^ 


+ m 2 (l -y - z) - ie 


(16) 


where A yz is given by 

Ayz = ym\ + zm\ - y(l - y - z)p 2 0 - z(l - z - y)p{ - yzpj . (17) 

Equation (16) is evaluated numerically and the real and imaginary parts are shown in Fig. [6] 


(left panel). If the widths of all intermediate particles are set to zero, the detailed structure of 
the amplitude becomes apparent. The imaginary part starts to grow rapidly from threshold 
yTftJh = uik + itik* and goes to infinity when yfs — E\. It exactly reproduces our result from 
section HI shown in Fig. [IJ The real part has a cusp at the threshold, then sharply drops 
below zero at a/s = E\ and becomes stable for higher values of y/s. 
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lm, Re 


Im, Re 




^V~s (GeV) 


Figure 6. (Left) Real (brown) and imaginary (orange) parts of A.^^ 0 (s) and M^ c 2 >7r j o (s), respec¬ 
tively. (Right) Real (brown) and imaginary (orange) parts of (s) when the finite width of 

K* is taken into account in VPP case. 

B. Realistic case: VPP intermediate particles 


For the realistic case of vector and pseudoscalar intermediate particles the expression for 
the matrix element is 

k ^-) (pi - h) u 

(18) 


[f( V PP) - 


d k\ 


£ 0n [ 9 


1 fM ' _ ^1 |~1 


(27t) 4 z (mf — k\ — ie){m\ — (p 0 — ki) 2 — ie)(m| — (k\ — pi) 2 — ie) 

We can apply the same procedure as in section IV A introducing Feynman parameters 
and performing a Wick rotation, and then calculate the resulting integral numerically. The 
details of the calculation are shown in appendix B 0 a Using the relation between 
and 


ai ->7T 


A ( V PP) 

^ai^Tr/o 


given in equation (14), the result is 

i r 1 r x -y i 


167T 2 

1 / 


d y 


d^ 


/ 0 A yz + m\(l — y — z) — ie 
'•i -y r 1 -y~ z 


iy .L iz L ix 


( y4po ■ Pi) + £pI _yy 


(19) 


\ (A yz + m\x — ie) 2 A yz + m\x — ie 

The real and imaginary parts of are compared to the scalar case in Fig. |b] (left). For 

both real and imaginary parts, the result for the VPP case is very similar to the scalar one, 


as was already shown for the imaginary part in section III 


C. Corrections to K* —> Kir vertex 


There are two additional corrections to be taken into account in order to arrive at a 
realistic estimate of the triangle amplitude: 
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1. Finite widths of intermediate particles. Until now we have assumed that the particles 
in the loop are stable (e —y 0). While this is reasonable for K , the width of K* is 
r x* = 0.05 GeV, 


2. P -wave tail suppression. In the VPP case the K* decays to Ktt in a P- wave, which 


is propagated to the / 0 7r final state, 


iy( v PP) 

Ui—>- 7r/o 


(Ai;t 0 ) \pA\ {\Pn\ is ai -)• foTT 


break up momentum). Therefore, the final / 0 and tt effectively are in a P- wave. This 
gives rise to an enhanced, unphysical tail in the signal intensity. 

To take into account the finite width of intermediate particles we substitute propagators 
of stable particles by resonance propagators. Technically this leads to substitutions m 2 —>■ 


m 2 — irrijTj. Including such a term for the K* propagator in equation (18) results in a 
smoother behavior of the amplitude, as shown in Fig. [6] (right panel). The singularity at 
y/s = Ei = 1.42 GeV is now limited and proportional to logT^*. 

For a two-body decay with orbital angular momentum L , the amplitude behaves like p L 
close to threshold due to centrifugal barrier. Far away from threshold, this is no longer 
correct because of the finiteness of the strong interaction. Accounting for the finiteness of 
interaction, however, is not unique. For a direct decay of a resonance phenomenological 
form factors are usually used, which come from a classical potential model. These could be 
Blatt-Weisskopf barrier factors |18] or exponential factors for finite meson-size corrections 
[ 19j . Another phenomenological approach is to introduce a simple left-hand singularity in 
the amplitude as a vertex form factor [20j. We demonstrate here that the latter approach, 
where a pole is introduced in the amplitude to account for the K* P -wave decay, gives a 


reasonable result. To do so, we include in equation (18) a factor 

F(h) = 


( 20 ) 


M 2 - k\ 

under the integral, where M is the position of the left-hand singularity and C is a constant 
normalized to the K * —> Ktt decay from mass shell, C = M 2 — m 2 Ki ,. Above the Ktt threshold 
this correction behaves like a .D-wave Blatt-Weisskopf factor ( F j^). So M corresponds to 
the size of K* 


1 + R 2 


1 + R 2 \p\ 2 « 1 + R 2 
(k\ - + m K ) 2 ) 


{kf - {m n + m K f)(kl - (m n - m K ) 2 ) _ 

Ak 2 

= M 2 = {m n + m K ) 2 — . (21) 




4 
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Figure 7. The intensity (left side) and the phase (right side) of the final amplitude A ai ^. n f 0 (s) 
including finite width of K* and P -wave tail suppression (see equation @) as a function of i/s 
for different values of the suppression parameter R. 


With this form of correction, our standard approach of Feynman parameters and Wick 
rotation can be used; the details of the calculation are again moved to appendix [Bj The final 
expression is 


^ / dy I dz x 

167r 2 y 


c 


"i -y 


(i -y- z) 


1 2 
4M 2 


(A yz + m\{ 1 - y - z)){A yz + M 2 ( 1 - y-z )) 

_(z 2 Pi + j/g(gi • P 2 )Kl - y-z) 2 _ 

(A^ + m\ (1 - y - z)) x (A^ + M 2 (l - j/ - ^)) x A yz 
1 . A yz + m?(l -y-z) 1 A yz + mf(l -y-z) 




m; 


A 




m? — M 2 


log 


( 22 ) 


A yz + M 2 (l - y - z))_ 

The final result including the final width of K* and the suppression of the P -wave tail at 
higher energies is plotted in Fig. [Tj The left panel shows the intensity, [A <Zl _ > . 7r / 0 1 2 (4p 2 /3) $ 2 , 
where $2 is the two-body phase space, for different values of the size parameter R = 
0.8 — 1.2 fm, and also without suppression (R — 0). The tail is indeed suppressed as ex¬ 
pected, almost independent of the exact value of R. The phase of the signal for R = 0 (no 
suppression) is shown in the right panel. 


Including the suppression factor in the integral of equation (18) artificially shifts the 


phase to lower values with respect to the case with no suppression. The phase motion, i.e. 
the relative difference as a function of the energy is not affected. Therefore we show only 
the phase for R — 0. 
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V. THE REACTION -k p —> a 1 (1260)p —> /o(980)7r p 


A. Cross section 


With a high-energy pion beam, as used in COMPASS and VES, the a] - (1260) is produced 
in a diffractive process proceeding via t-channel Pomeron exchange between the beam 7 r~ 
and the target proton, as shown in Fig. [ 8 j 


7T-' Pi 



Figure 8 . Diagram for diffractive production of the ai(1260) by scattering of a high-energy tt~ 
off a proton target, which remains intact. The tt~tt + tt~ final state observed in the experiment is 
modeled by the decay of the into a charged pion and a neutral isobar R (p° or /o(980)), which 
subsequently decays into two charged pions. The isobar can be formed either by or by ir vr 3 ". 


In order to estimate the intensity of the signal expected in the /o( 980 ) 7 t“ channel we 
calculate its intensity and phase difference compared to the dominant a]” (1260) —» p° 7 r~ 
decay, assuming that the signal in / 0 7 r _ is entirely due to the triangle singularity in the 
decay ab(1260) —» / 0 ( 980 ) 7 t _ . We denote the invariant mass squared of the a\ by s, and 
the isobar invariant mass squared by Si 2 or s 2 3 , respectively. Factorizing out the production 
cross section cr pro d(s) of the ai(1260), which is independent of the final state, the differential 
cross section for the full process n~p — y a)”(1260)p — > Rn~p — y 7i~7i +/ n'~p can be written as 


do" O"p roc i(s) 

ds 47 t 


dsi 2 2 m ai r ai (s, S 12 ) 2 mijr R ^. 7 r + 7 r -(si 2 ) 

2 vr (m 2 ai - s ) 2 + ( s , s 12 ) (m 2 R - s l2 ) 2 + m 2 R T 2 R (s 12 ) 


+{12 o 23} + 


d$ 37r Interf(si 2 , s 23 ) 


(23) 


where T a ^b c is the partial width for the two-body decay a —» be, m aii R and T aijR are the pole 
masses and mass-dependent full widths of a± and R , respectively, and d<h 3ff is the 37 t phase 
space. The first two terms constitute the contributions with the isobar in the and 

in the 7 h} subsystem, respectively. The third term is the contribution of the interference 
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between two processes. The latter is found to be very small for oy —> p°tt~ (less then 2% for 
oy decay from its mass shell) as well as for ay —> / 0 7r _ , so we will disregard that contribution 
and use the following equation for the cross section: 


d(T ^ O"prod( , 5) 
ds 2tt 


ds 


12 


2m ai r o -^ R7r _(s,s 12 ) 




2vr (m 2 ai - s) 2 + m2 T2 (s, s 12 ) {m 2 R - Si 2 ) 2 + m 2 R T 2 R (s 12 ) 


(24) 


The partial widths of ay and isobar decays are calculated by averaging the expressions 
for the square of the matrix elements of the corresponding hadronic vertices over the initial 
states and summing over the final ones and multiplying by the corresponding phase space. 
For the decay ay (1260) —* p°n~ (axial vector to vector and pseudoscalar, AVP), we have 


T ar ^ p07r -(s,m 2 ) = 


2 m r 9a iP°' K 


°CL1 


i + 


3 m 2 


x 


1 2\p p \ 
8tt ’ 


(25) 


while for oy (1260) —> / 0 vr (axial vector to pseudoscalar and scalar, APS), we get 


a l ->fon 


- (s ’ ro2)= 2^ 9 b..- (s ’ ro2) 


4 IP/o | 2 w 1 2 \Pfo\ 


X 


87T y/s 


(26) 


where \p p /f 0 (s, m 2 )\ = X 1 ^ 2 (s,m 2 ,m 2 )/(2y/s) is the break-up momentum for the two-body 
decay of a particle with mass y/s to particles with masses m = yfsvi and m,. The coupling 


of oy to foTT in equation (26) is given by 

g 2 - f _(s,m 2 ) = 

Jon v ’ ' 


. ( v pp) 

^ai^fow 


2 / ( K*K+c.c ,) N 2 
9 a i fon~ 


(27) 


where ^ PP Jf Q7T is the triangle amplitude calculated in equation (22) and / +c - c -'> j s an 


a i- 


effective coupling which includes the couplings of the individua 
diagram, taking into account both isospin channels. 

The expressions for the isobar decays are 

1 4|p 7r | 2 


«1 fo™ 


vertices in the triangle 


bpO—»7r+7T- (gn ) 


2m, 


2 

9pO n + ir~ 


x 


1 21 p n 
87T 771 


and 


1 2 

i / 2\ L ^ 2 

/o ^ +7r -(m ) - ^^ 3 


1 2\p n \ 2 _ 

9f °™ X 8tt m ~ 3^ 0 ™W ’ 


(28) 


(29) 


where the dimensionless coupling g/ ovr7r = fl , / 07r has been introduced, and |p„-| = 
A 1 A(tt 7, 2 , m 2 , m 2 )/(2m) is the break-up momentum of the isobar with mass m to two pions. 
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B. Evaluation of the couplings 


To evaluate the magnitude of oq —>• / 0 7r decay with respect to cp —* pn S -wave we take 
into account two possible isospin configurations of intermediate states (. K*°K~K + ) and 
(. K*~K°K °) and evaluate the corresponding couplings. The table gives the couplings and 
Clebsch-Gordan coefficients for each vertex inside the loop of the two isospin configurations. 


Vertex 

K*°K~K + 

K*~K°K° 

a i 

g ai K*K 

g ai K*K 

n~ 

\/2/3 9k*Ktt 

\/ 2 /3 9k*Ktt 

f° 

V 71 / 2 9 foKii 

y/ 1 / 2 9 foKii 


Since the isospin structure of both configurations is identical, the two diagrams add up. 
Disregarding the mass difference between charged and neutral kaons, the contributions of 
both are the same, so the effective coupling of process can be written as 

2 


9 


(. K*K+c.c .) 
a i foir- 


~ ^y^9aiK*K 9K*K-n 3f 0 KK 


(30) 


We first consider a\ decays. The resonance is rather wide, so the energy dependence of 
the width should be taken into account. The best knowledge about a\ decay channels and 
branching ratio comes from hadronic r decay measurements pa eh ezi . For simplicity we 
consider only the main contribution to the energy dependence of the width, which comes 
from ai —> pit S- wave. 

\p(s)\ m ai 


r ui (s) = r ai (m ai )- 


aiV ~' ^ 1 v "" i; Ip«)I 


s = 


A 1,/2 (s, mV mV 


p’ 


(31) 


2^ 

We use the measured branching ratios, Br(Gq —» pn, S' —wave) ~ 60%, Br(ai —» K*K + c.c., 
S — wave) ~ 2.2% to extract the ratio of couplings. The ratio g aiP n/g ai K*K is calculated 


with the help of equation (24), where the energy dependence of the production mechanism 
is disregarded. For «i —> pi:, a size correction form factor exp(— K 2 \p p | 2 ), where p p the is 
break up momentum, is applied [TS]. For R > 0.5 GeV -1 , the convergence of the integral 


over ds in equation (24) is achieved for an upper limit of the integration < 5GeV, while 
for R = 0 a higher limit is required. Varying R between 0 and 5GeV _1 , and including the 
uncertainty due to the slow convergence for R = 0, the resulting ratio of the couplings is 


laipiv 


2 gV 0 

CL i P U 7T 


9. 


Oi K*K 


gt- 


K*°K- 


+ 9„- T/*- 


6-10 


(32) 


K*~K° 
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For the evaluation of relative strength of f 0 n signal in section VC, we use 9 2 iP -k/ 9 2 aiK *R = 


The width of K* is measured precisely and the branching to Kn P -wave is 100% [T], The 
corresponding coupling can thus be extracted from 


IV = 


1 


2m jo 


2 4|fc| 2 1 2|/c| 

* o-~ > 

3 on rrif<c* 


\k\ = A 1/2 (s, m 2 K *,mD/(2m K * 


(33) 


to be — 31.2. 

The parameterization of the /o is not trivial, since both decay channels (/o —> 2n , /o —>■ 
2 K) need to be taken into account. We make use of the Flatte parameterization (23] of the 


/o propagator and the decay width in equation (24), 

2m /o VrlVI 


- S 12 - im fQ (g nn \p n \ + g KR \pi< 


, \Pn/ K {s 12 ) \ = ^s l2 -Aml /K 


(34) 


The measurements of the branching ratio Br(/o —> KK)/ (Br(/o —> nn)) and the ratio of 
couplings extracted therefrom, Rk/-k = 9j 0 kR/^ fom, = 9f Q K k/ 9fo-K-n ~ 4, are rather consistent 
with each other [241 25] . But the absolute values of the couplings are not known very well. 
For our estimation of the branching we use g/ 0 ^ = 2.3 GeV [25] , so 9f Q7rR = 0.21, gf 0 KK ~ 0.8. 


C. Evaluation of the branching ratio 


The cross section calculated with the help of equation (24) for cq (1260) —> fan —>■ 
n~n + n~ is compared to the one for the dominant channel < 2 ^( 1260 ) —)■ p°n~ —> n~n + n~ in 
Fig. |9j Here, the peak of cq —>■ pn has been normalized to 1, the f 0 n channel is shown in 
relative scale. Under the assumptions detailed in the previous subsection for the couplings, 
the peak-to-peak ratio is 1 : 100, in very good agreement with the experimental result. 

Let us discuss how reliable our estimation is and which factors could affect the magnitude 
and the shape of fan peak. First of all, we assumed that the origin of both K*K, which 
are rescattered to f 0 n, and pn is decay of the cq(1260) resonance. In a hadron fixed-target 
experiment like COMPASS or VES, however, there may be other processes which contribute 
to the same final state, e.g. non-resonant Deck-like processes IZH- We expect a rather large 
contribution of Deck-like background to pn S -wave signal EU as well as to K*K channel 
J. Taking into account these processes, the resulting ratio of couplings could be different 


from equation (32). 
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CTa, (a.U.) 



Figure 9. Cross sections for ai(1260) resonance as a function of invariant mass y/s in arbitrary 
units. (Left) comparison of dominant a—y p 0/ K~ S -wave decay and a± —y P -wave channel 

due to rescattering of kaons. (Right) Pseudo-resonant shape of a^ > /qvt^ in linear scale. 


Secondly, pn —y f 0 n rescattering plays an important role. Using the same method as in 


section IV B one can show that the triangle diagram a\ —> p7r —» f 0 n gives a rather flat 
amplitude with a constant phase (magnitude ~ 4% of peak intensity of K*K + c.c. —> / 0 7r). 
It interferes with the signal from K*K+c.c. and changes its intensity and phase. Taking into 
account this contribution is in principle possible, but requires the knowledge of the relative 
sign between a\ —> pix and a\ — > K*K, i.e. the relative sign of ss in ai, which is unknown. 

The third uncertainty comes from the shape of the f 0 and the corresponding coupling 
constants. We found that the shape of our signal is stable for different values of gf on7r and 
Rk/tt■ The relative intensity, however, is proportional to n, which could therefore 

easily change by a factor of two depending on the input values. 
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VI. CONCLUSIONS 

Even after many years of intense studies, both experimentally and theoretically, the 
excitation spectrum of hadrons is still not understood. This is especially true in the region 
of charm and bottom quarks, but also the light-quark sector sometimes bears surprises. 

Present-day experiments are collecting extremely large event samples, with allow them 
to perform analyses with very small statistical uncertainties and permit them to find small 
signals which were not observable before. Recently, the COMPASS experiment has reported 
the observation of a resonance-like signal with axial-vector quantum numbers J FC = 1 ++ in 
a completely unexpected mass region only about 0.2 GeV above the ground state ai(1260), 
decaying to /o( 980 ) 7 t. 

In this paper we show that a resonance-like signal with a maximum intensity at 1.4 GeV, 
compatible with the experimental result, can be generated dynamically via a triangle singu¬ 
larity in the decay of the ground state cq(1260) to K*K + c.c and the subsequent rescat¬ 
tering of the K from K* decay to form / 0 (980). This process also generates a rather sharp 
phase motion, which is not locked to the phase of the wide Oi(1260). The singularity appears 
in a kinematic region where the intermediate particles are collinear and on mass shell. The 
structure of the amplitude is investigated in two ways: first, the imaginary part is calculated 
using cutting rules; second, the full amplitude is evaluated using Feynman rules in order to 
obtain the imaginary and real parts. Both approaches are performed for the hypothetical 
case of scalar intermediate particles and for the realistic case of vector and pseudoscalar 
intermediate particles. It is shown that both cases give very similar results. For the final 
result of the triangle amplitude, we also include the finite width of the K* and a phenomeno¬ 
logical factor to suppress the tail clue to the P -wave decay of the K* at high energies. The 
inclusion of this factor, however, is not unique, and also spoils the beauty of the solution 
somewhat, because it shifts the phase of the amplitude which is questionable and needs 
future investigation. The treatment of decays with higher angular momenta inside loops in 
terms of analytical solutions certainly needs more attention in the future. 

We then estimate the magnitude of the signal expected in the / 0 (980)7T channel due to the 
triangle singularity compared to the dominant decay of the ai(1260) to p7r, also observed 
in the experiment. Our result gives a relative peak intensity of 1% for the foil channel, 
with a rather large uncertainty which is due partly to poorly known couplings and partly 
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to other rescattering processes like a± —>■ pn —> foir. The last process does not produce a 
singularity in the kinematically allowed region, but the corresponding amplitude interferes 
with the K*K amplitude and modifies its intensity and phase. In addition, we only consider 
the genuine ai(1260) resonance as source for the triangle diagram, while it is known that 
in the reaction ir~p —y 3np there is a rather large contribution to the intensity in the pi r 
channel from non-resonant processes like the Deck effect, which may influence the relative 
branching ratio. 

The dynamical interpretation of the ai(1420) presented in this paper captures the main 
effect and probably accounts for a large fraction of the signal observed by COMPASS and 
VES. As a next step, one may fit our amplitude to the data and compare to the Breit-Wigner 
fit, and eventually extract better values for the coupling constants. The data sample on the 
ai from r-decays should also be large enough to observe the f 0 n peak if the data is fitted 
without phase locking of f 0 n with a\. In general, the large data samples available nowadays 
both for light and heavy hadrons allow us to revisit effects which were already discussed more 
than 30 years ago, but were almost forgotten since then because data were too scarce to test 
them. These may play an important role in our understanding of the hadron spectrum. 
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Appendix A: Parameterization of vertices and Feynman rules. 


First we mention the approach we use to parameterize the vertices for interactions of 
particles. From symmetry considerations the Lorentz structure for vertices is the following 
(we use S = scalar, P = pseudoscalar, V = vector, A = axial vector): 


£o, Po 



= 9vpp - kzY 


Po 


S 

►—< 


P 

* 



= 9spp 


Figure 10. Parameterization of VPP vertex (P-wave) and SPP vertex (P-wave) 



9AVP £q £ i p 


Figure 11. Parameterization of AVP vertex (S- wave) 


ki, mi _ 1 

P m^—k^—ie 


h,mi _ g^-k^k\/kl 
V nij—kf—k 


Figure 12. Propagators for pseudoscalar and vector particles 


Appendix B: Calculation of integrals 


In this section we calculate three integrals: 

d k\ 


V« = 


(27r) 4 f A 1 A 2 A 3 ’ 


Va = 


d kt £ ^ (■<" - ^r) (pi - ^'3) 1 


V* = 


J (2tt) 4 i 

dki £ om (V" -n^)(Pi-h)vX J IP 


AiA 2 A 3 

fcffch 


j (2vr) 4 i AiA 2 A 3 

where Ai = mf — A; 2 — ie, A 2 = — (p 0 — &i) 2 — *e, A 3 = — (&i — pi) 2 — ie. 


(Bl) 

(B2) 

(B3) 
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a. First integral 


For the calculation of V 3 , Feynman parameters (x,y,z) are introduced to rewrite the 
integral: 




14 = 



dx d y cb 2! S(x + y + z — 1) 


d A k 


(B4) 


'0 Jo Jo J (27t ) 4 iD 3 

where D = x(m\ — k\ — ie) + y(m\ — (p 0 — k \) 2 — ie) + z(m\ — (£q —pi ) 2 — ie). The quadratic 

form D{k\ ) can be reduced to diagonal form collecting terms with k\ and extracting the full 

square. The condition x + y + z = lis used. 


D = —(A', -yp 0 
■2 


zpi ) 2 + A - ie, 


A = xmf + A yz , A yz = yrrq + zm\ - y(l - y - z)p „ ~ z(l - z - 


y^p 2 m 


After shifting the variable of integration k\ —* / = Aq 
V 3 — [ [ [ dxdydz 2! 5(x + y + z — 1) 


ypo — -spi we have 
r ' d 4 / 1 


(B5) 

(B6) 

(B7) 


jo jo jo J ( 2 ir) 4 i (-; 2 + A-K) 3 

For the integration over 2 0 , notice that the denominator has poles when Iq — r + A 2 — te. 
The positions of the poles are functions of the external invariants p\, p 2 , p| and the Feynman 
parameters. The basic idea, which we use is aimed to calculate the integral in the region 
where pi < 0, p\ < 0, p\ < 0, i.e. A > 0 for all value of x,y,z. In that region we can 
rotate the contour of integration over Iq anticlockwise (Wick rotation) and integrate along 
the imaginary axis. We make use of the transfer of the integration variable l to Euler space 
with integration variable Ie [} 2 = —1%), where the integration is much simpler. One has: 

d H e 1 




w = 



dx d y dz 2! 5(x + y + z — 1) 


1 0 Jo 


(27t) 4 (/| + A - ie ) 3 ' 


The next step is integration over Ie- As a result we have: 

1 r 1 r l ~y 

iz 


A„* + m 2 (l - y- z) - it 


(B8) 


(B9) 


Equation (B9) is simple enough for numerical integration. 


b. Second integral 


Let us consider integral V 4 , equation (B2). The numerator can be simplified as 

i-i'i- 1 ' 

/v-j^ /v-j^ 


£ 0fi 9 


/j,u 


k 2 


(p\ - h) u = e 0fl g 


fiv 


k i k i\ , n Of \ 1 O ( £ o • h)( k i - Pi) 

( 2 p 1 -k 1 ) u = 2 [£q - pi) + 2- 


k l 


-kl 


(BIO) 
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One can notice that kf in the numerator has the same form as Aq = — kf with mass 


nig = 0. The integral in equation (B2) is equal to 


- 0 4 = (£ 0 -Pi) 


d 4 Aq 


+ 


d 4 fci (e 0 • fci)(fci ygi) 

(2v t)H A 0 AiA 2 A 3 


(Bll) 


(27r) 4 7 AiA 2 A 3 

The first integral in equation ( Blf[ ) is equal to equation (B9). For the second one we 
introduce four Feynman parameters: 

/' d 4 Aq (e 0 • fci)(fci-pi) 


(27r) 4 i AoAiA 2 A 3 


dt 


dx d y dz3\ 5(t + x + y + z — 1) 


d 4 ki (e 0 • fci)(fci vgi) 

(27t) 4 i D 4 


, (B12) 


where for D 4 with condition xq + X\ + x 2 + X 3 = 1 we have the same expression as equa¬ 


tion (B5). So the same shift of the integration variable k\ is used, i.e. Aq —> l = ki~yp 0 — zpi. 
The expression in the numerator can be written as 

(e<i-ki)(h-pi) = [coPi] + l p-[( z Pi+y(Pi-P 2 ))eo+z(s„-p 1 )p^]+z(e 0 -pi)( 2 pl+y(pi-p 2 )). 

(B13) 

After Wick rotation and the integration over angular variables dfl 4 , the term proportional 
to gives zero and —> —g yv l 2 E /A. So one arrives at 


d 4 fci (g 0 • fci)(fci -pi) 

(27r) 4 i AoAxAaAg 

1 f dH K 


= / d t dx 


x (e 0 • Pi) 


J 2 

l E 


4 J (2tt) 4 (ll + A - ieY 


dy dz3\ 5(t + x + y + z — l)x 
Jo Jo 

+ z(zpl + y(pi -p 2 )) f dlE 


(27t) 4 (1% + A - ie) A _ 

(B14) 


All integrals converge. The integration over dt is removed by a delta-function. 


Im - ( g ° 'P^ 

2 4 16tt 2 


’ 1 - 2 / 


L/o 


dy 


dy 


" 1 - 2 / 


d^ 


A yz + m\{ 1 -y - z) -ie 


-l-y-z 


d^ 


dx 


( yzjpo-Pi) + z 2 p\ _ 1/4 


\ (A yz + m\x — ie ) 2 A yz + mix — ie 


(B15) 


A yz is given by equation (B 6 ). 


c. Third integral 


The calculation of V 5 , equation (B3) proceeds similarly to V 4 . The difference is that we 


have four poles instead of three in the denominator of equation (B3). First, one can simplify 


























the numerator as in equation (BIO): 
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1 f d 4 k 1 1 r d 4 Aq (e 0 ■ fci)(fci -pi) 

2C 5 " l£ °' PlJ j (27r) 4 iA 1 A 2 A 3 A 4 J (2n)H A 0 A 1 A 2 A 3 A 4 ’ 


(B16) 


where A 4 = M 2 — k\. 

We introduce four and five Feynman parameters for the integrals, respectively. The 
expression for the denominators are .Dg and D- 5 : 

D 5 = x(m\ - k\) + ?/(m 2 - (p 0 - £q) 2 ) + z(m\ - (fc 4 - p 4 ) 2 ) + u(M 2 - A; 2 ), (B17) 


= t(-A; 2 ) + x(m\ - k\) + - (po ~ h) 2 ) + z(m\ - (k x - p 4 ) 2 ) + u(M 2 - k\). (B18) 


Then we perform the same calculation as in section B 0 b, and the result is 


J_ r = ( g o ’Pi) 

2 C 5 16vr 2 „ 

rl—x—y—z 


d y 


" 1 - 2 / 


"1 -y-z 


dx 


dx 


d u 


yz(p 0 -pi) +z 2 pl 


(A yz + m\x + M 2 (l — x — y — z) — ie) 2 
1/4 


(A yz + m 2 x + M 2 u — ie) 3 (A yz + m 2 x + M 2 u — ie) 2 


, (B19) 


where A yz is given by equation (B6). To make the expression simpler and convenient for a 
numerical evaluation, we carry out the integration over dx and d y explicitly: 


2 C V5 16tt 2 


dy 


-l-y 


dz X 


(l-y- z) 


(■ A yz + ml (1 -y-z)- ie)(A yz + M 2 (l - y - z) - ie) 
_ (z 2 p\ + yz(pi • p 2 ))(l - y-z) 2 _ 

(A yz - m'i( \ - y - z) - ie){A yz + M 2 (l - y - z) - ie)(A yz - ie) 


1 2 
4 M 2 


1 A yz + m\(l - y- z) - ie 

Z2 lo s - 


m 


l A^ z ^ 


771T 


M 2 


log 


A yz + ml(l - y - z) - ie 
A yz + M 2 (l -y - z) -ie)_ 


(B20) 
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